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Abstract 

The two point integrals contributing to the self energy of a particle in a 
three dimensional quantum field theory are calculated to two loop order in 
perturbation theory as well as the vacuum ones contributing to the effective 
potential to three loop order. For almost every integral an expression in 
terms of elementary and dilogarithm functions is obtained. For two integrals, 
the master integral and the Mercedes integral, a one dimensional integral 
representation is obtained with an integrand consisting only of elementary 
functions. The results are applied to a scalar A0 4 theory. 
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1 Introduction 



When physical phenomena are described by quantum field theories, all the 
observable quantities are expressed in terms of functional integrals. Since 
these integrals can be evaluated exactly only in very special cases, one has 
to use lattice simulations or analytical approximation methods. The most 
common method is perturbation theory. One expands the desired quantity 
as a series of integrals represented by Feynman diagrams. 

Although the field theories of particle physics are four dimensional, the 
importance of three dimensional theories has grown recently. The main rea- 
son is the method of dimensional reduction of a four dimensional finite tem- 
perature field theory to a three dimensional zero temperature effective theory 
[1-6]. This technique has been applied to the electroweak phase transition of 
the early universe [7-9]. Three dimensional field theories are also important 
in the theory of critical phenomena. 

The nature of the phase transition is an important question in both main 
applications of three dimensional field theories. Thus the effective potential 
which gives the true ground state of the system has an essential significance. 
Unfortunately, in dimensionally reduced effective theory, perturbative calcu- 
lations are applicable neither in the symmetric phase nor in the immediate 
vicinity of the phase transition. However, the perturbative results obtained 
deep in the broken phase can give new insight into the problem and the 
lattice results [§, §]. 

At one loop level the perturbative calculations are fairly easy. When 
one needs higher corrections, the integrals get more complicated. In four 
dimensions these integrals have been studied recently by many authors [10- 
19] but a three dimensional discussion has been missing. 

The purpose of this paper is to evaluate all the integrals needed for the 
self energy of a particle in a three dimensional scalar theory to two loop order 
and the ones needed for the effective potential to three loop order. As shown 



by Weiglein et al. [20"], the self energy of a particle in a gauge field theory 
can be expressed in terms of these scalar integrals. Most of the integrals are 
calculated in a straightforward way, but with two integrals a different route 
must be chosen. 

The paper is organized as follows. In Sect. ^| the simpler integrals are 
calculated explicitly. Sect. |3| is devoted to the two more intricate integrals. 
In Sect. £| the integrals are applied to the calculation of the self energy and 
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Figure 1: Topologies of one and two loop vacuum diagrams and one loop two 
point diagram 

the effective potential of scalar A0 4 theory. 

2 Evaluation of scalar integrals 
2.1 Classification of integrals 

In order to calculate the self energy of a specific system one needs to consider 
all possible one particle irreducible two point diagrams. Let us assume that 
the Lagrangian consists only of terms at most quartic in the fields. Then the 
diagrams can be composed from three and four leg vertices using Feynman 
rules. 

Depending on the number of fields and the specific form of the Lagrangian 
the number of possible diagrams may vary. However, there is only a restricted 
set of different topologies these diagrams may have. At one loop level there 
are only two possible topologies. These two are the diagrams a and (3 shown 
in Fig. |l|. At two loop level the number of different topologies is eight. These 
are given in Fig. |2|. 

For the effective potential one needs the vacuum diagrams |T2TJ. At one 
loop level there is of course only one possible topology, and two at the two 
loop level. At the three loop level the number is six. These diagrams are 
shown in Fig. [3[ 

Some of the integrals corresponding to these diagrams do not converge, 
but when they do, there is a relation between vacuum and two point integrals. 

Lemma 1 Let /2G0) be a two point integral with an external momentum of 
p and F (m) the vacuum diagram obtained by connecting the outer legs of f^ 
with a particle of mass m. Then the following relation holds: 
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Figure 3: Topologies of three loop vacuum diagrams 

The lemma can easily be proven by taking a Fourier transform of both 
sides of Eq. (0). 

One should also notice that in case there are more than one propagator 
with the same momentum, they can be separated to two integrals by partial 
fractioning or be written as a derivative: 

1 

(p 2 + m 2 )(p 2 + 777J) 
1 

(p 2 + m 2 ) 2 

Some of the integrals also factorize into separate parts. Using these rela- 
tions the number of integrals to be calculated can be reduced to six: a, /3, 
g, A, C and E. Of these, diagrams a and (3 are easily evaluated. Diagram 7 
is only a special case of g. Diagrams a, b and f are related to diagrams A, B 
and C, respectively, by Lemma [I]. Finally, all the diagrams c, d, e, h, D and 
F factorize to products of simpler diagrams. The results for all the integrals 
are given in explicit form in appendix [A]. 
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2.2 Diagram g 

Consider now the diagram g called the sunset. Diagram 7 is a special case 
of this with vanishing external momentum. Let the masses of the particles 
be mi, m 2 and m 3 . In this integral the coordinate space method P^ , |23"| ] will 
be used. In 3 — 2e dimensions the Fourier transform of the propagator is 

Vi{x) = (Vr— L- i^-Y' 6 K,(m lX ). (4) 



(2vr)^- £ V x 



2 



Here K u (x) is the modified Bessel function. After the transform the integral 
reads 

I g (k;m u m 2 ,m 3 ) = /i 2e / d^Re^H^R). (5) 

i 

The divergence occurs only on the limit R — > 0. Thus the integration can be 
split at R = r: 

X g {k) = (^fj ' £ dRR*-<J^{kR)Y[V i (R) 

A-jr POO 

+— / dRRsm(kR)Y[Vi(R) 

k Jr i 

= T { ;\k)+if\ki (6) 

where J v (x) is the Bessel function and p, is the MS scale parameter fx 2 = 

Since converges, the usual three dimensional Fourier transform of the 
propagator can be used. In X^ a \ that is when R < r, one can approximate 
the Bessel functions by the lowest order terms in their Laurent series: 

Jj_(W2) = f ^—^(lkRy~ e + 0((kR)i), (7) 

vtm ^Vy r(|- £ ) 1 l+2e 

-(eV) e r( r ( Lt ) e) ^ 2e + o(fi)- (8) 

The error vanishes as the limit r — ► is taken. For the present integral only 
the 0{R~ 1 ) term of Vi is needed. Now one is left with only a straightforward 
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task of integrating over powers of R. All the terms that are singular at the 
lower limit R = are treated by analytical continuation to sufficiently great 
values of e so that they vanish. 
The result is 

^W = ^(^ + log^ + i + 7 ). 0) 

T*p is a normal three dimensional integral and it can be calculated using 
normal methods of multidimensional integration: 

Z?{k) = -Lfl- mi+T? ; 2 + m3 arctan k 



(An) 2 V k mi + m 2 + m 3 

-7 - logr - - log ((mi + m 2 + m 3 ) 2 + k 2 ^j\ . (10) 



2 

Now one can write down the result: 

1 / 1 mi + m 2 + m 3 k 
X ^ k ) = 777^7 77 1 arctan ' 



(Air) 2 \4e k m x + m 2 + m 3 

1 ^ 3\ 

+ 2 ° g (mi + m 2 + m 3 ) 2 + A; 2 + 2/ ' ( j 

2.3 Reducible diagrams 

The integral g together with the well known results for the integrals a and (3 
make it possible to calculate all the reducible integrals using Lemma [1] and 
Eq. (fj). As a simple example, consider the diagram c: 

d 3 ~ 2£ k d 3 - 2£ q 
(2tt) 3 - 2£ (2tt) 3 - 2£ 

' (12) 



J, 



ik 2 + m 2 )(q 2 + m 2 )(fc 2 + m 2 )((p — k) 2 + m\) ' 
The g-integration is nothing but I a . One then uses Eq. (g) to obtain 
m 2 1 r d 3 ~ 2e k ( 1 1 \ 1 



J c 



47r m 2 — mf 7 (27r) 3 ~ 2e \ k 2 + m 2 k 2 + m 2 J (p — k) 2 + m 2 
m 2 1 



47T m 2 — m? 



(J/3 (mi, m 4 ) - J^(m 3) m 4 )) 



arctan arctan ) . (13) 



(47r) 2 p(m 2 — m 2 ) V mi + m 4 m 3 + m 4 
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In a similar way the diagram D can be factorized: 

1 D = TJm^)—^- — r (lJmi,m 2 ,m 3 ) - lJm 1 ,m 2 ,m i )) 

m 5 m 1 + m 2 + m i 

log ■ ■ . (14) 



(47r) 3 (???J — m|) m\ + m 2 + m 3 

The results of the integrals are convergent, but the integrals themselves 
are actually not, since they both contain the simple loop a, which is divergent. 
In dimensional regularization this divergence vanishes, but Lemma [I] not 
proven to hold with dimensional regularization. However, this loop is in 
both cases factorized as a separate integral and the remaining other factor is 
convergent. Therefore Lemma [l] should hold for these diagrams: 

l c (p;m 1 ,m 2 ,m 3 ,m i ) = 
2-ni 

— (I D (ip, m 4 , mi, m 3 , m 2 ) - X D (-ip, m 4 , m x , m 3 , m 2 )) . (15) 
V 

Substitute the previous result to the right hand side to obtain 

X c = - — r^— -. — | 7y- ( arctan arctan ), (16) 

(47r)^o(mf — m 3 ) \ mi + m^ m^ + rrii/ 



which is the same as Eq. (|i~3|). 



2.4 Diagram E 

Since the integrals E and g do not converge, Lemma [3] does not hold. There- 
fore one must start the calculation of E from the beginning. The coordinate 
space integral is 

r 4 

J B (mi, m 2 , m 3 , m 4 ) = fi 2e J d 3 ' 2e R ]J V^R). (17) 



n=l 



Just like before, this is separated to two parts, only one of which is di- 
vergent: 



4 (4vr 



L_^^i_ e -E. m ,H (18) 
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In the expansion of V{ one now has to take into account also the term of 
order 0(1). Then the first integral reads 



X 



(a) 



i / e yr £ r(|) 



(4tt) 



T(l-e) Jo 



dRR 



2-2e 



r( f -ej 



'eV' 



V r(|) 
^ /rfi- £ )\ 3 



R 



(eV) £ 



V 4 ; V r(|) 

This simple integral gives 



r(-| + gi 
r(-|) 



-3+6e 



E 



m 



l-2e 



(19) 



(a) 



1 



(47T) 



11 /l r 2 u 3 ' 

-E m M-+ 4 + 4 7 + 21og-^ 

r 4 ^ \ £ 2rrij 



(20) 



X^ is also easily evaluated and one can write the result for the whole 
integral: 



X K = t-^t jTm; -— + 2 + - log — + log — ^ 
\ 4e 2 // ii 



(An) 3 



(21) 



i=l 



2.5 Diagram C 

Consider now the diagram C. This integral converges and Fourier transform 
can be used directly in three dimensional space. The masses are as shown in 

Fig. a. 



x 



c 



1 



P ,q,k (p 2 + mf)((/c — p) 2 + mDdk — q) 2 + m§) 
1 



(q 2 + ml)(k 2 + m 2 ) 
Taking Fourier transform of this gives 



X 



c 



(4tt) 



d 3 x±- 



-(mi +7712)11 

x\ 



d 3 X2- 



-m 5 \xi+X2\ — {m3 +m4)X2 

\x\ + x 2 \x\ 



(22) 



(23) 
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Concentrate now on the latter integration 

g— m 5 \xi+X2\ — (iri3+m4,)x2 



d x 2 - 



\X 1 + X 2 \X2 



2vr 

+ (e m5Xl 
2tt 



c i dx 



x 



— e 



-(m3+rn4+m 5 )x\ 



poo fjqr 

e -m 5 x 1 \ / ~ c -(m 3 +m4,+m 5 )x 



m 5 xi 



log 



xi X 
1713 + 7774 + 7775 
7773 + 7774 — 7775 



+ Ei ((7773 + 7774 — 7775)^1 



-e msa:i Ei((7773 _|_ m4 _|_ m5 ) Xl ) 



(24) 



where Ei(x) is the exponential integral. 

Substituting this into Eq. and integrating gives the result 



c 



lim 



1 



7773 + 7774 + 7775 

log — — — (-7 - l0g(777i + 777 2 + m 5 )r) 



(47r) 3 27775 L 7773 + 7774 — 7775 
-- (((2) + (7 + log(777 3 + 777 4 - 777 5 )r) 2 ) - Li 2 

+ i (((2) + (7 + log(m 3 + 7774 + 777 5 )r) 2 ) + Li 2 



777! + 777 2 + 777 5 
7773 + 7774 — 7775 
TUX + 7772 — W5 



7773 + 7774 + 7775 / . 



(47r) 3 2m E 



. 7773 + 7774 + 7775, v /(777 3 + m 4 ) 2 -m| 

log ; log 



7773 + 7774 — 777 5 



+Li 2 (- mi+m2 ~ m5 )-Li2 

TTls + 7774 + 7775 



777i + 7772 + ^5 
777i + 7772 + ™5 
7773 + 7774 — 7775 



(25) 



Li2(x) is the Euler dilogarithm function and ((x) is the Riemann zeta func- 
tion. 



3 Mercedes and the master integral 
3.1 Mercedes integral 

The last diagram to be considered is A. Many different techniques have been 
proposed for the calculation [10-19], most often in four dimensions. The 
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Figure 4: The momenta and the masses of the Mercedes diagram 



coordinate space method of the preceeding subsection does not work well 
even though the integral converges and three dimensional transform can be 
used. The transformed integral is namely as complicated as the original one. 
Other approachs are based on some series expansion of the integral but then 
the resummation of the series is usually not possible. Many other methods, 
like Mellin-Barnes transformation and Gegenbauer polynomials, work well 
only in four dimensional space. In four dimensions one has been able to 
reduce the general combination of masses to a one- dimensional integral. 



The approach to be used here is similar to that of Kotikov p4j . A differ- 
ential equation which the integral must satisfy is constructed and then this 
equation is solved. For propagators a shorthand notation will be used: 

A ; = T7 — 2" ( 26 ) 

With the masses and the momenta defined as in Fig. ||, the integral reads 

l A (m 1 ,m2,m 3 ,m4,m5,m 6 ) = J AjA^Aj^Aj^AjAg 

A 123456 . (27) 

Here integration over all momenta is assumed. 

The first step is to integrate this by parts. The boundary terms vanish. 

1a = 



3 J dqi 
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-~ J (-2g 2 Aj - 2q ■ (q- k)A 2 q _ k - 2q- (q- p)A 3 p _ q ) A 



123456, 



(.28) 



Using conservation of momentum and then simplifying one gets the result 



1a 



J (2m\A\ - (A\)-'Al_ k + (A 6 k y 1 A q _ k + (ml + m\ - m 2 )A 2 _ fc 

-(K)~ lA l- g + (K^K-o + (ml + ml- ml)Al q ) A 123456 
2m 2 i^ 2 - J\ + J 2 + (ml + m\- ml)K 2 

- J 3 2 + Jf + (ml + ml - ml)K 3 . (29) 



Here K t is Ta with the propagator of particle i squared and J/ is a diagram 
of type If? which has been constructed by removing the propagator of j from 
Ia with propagator i squared. Then by Eq. (H) one has 



Ki 







2m,; dm; 



Z A {mii m 2 , m 3 , m 4 , m 5 , m 6 ) 



(30) 



J/ can be obtained from the following expression by substituting the corre- 
sponding masses 

2^ (mi, m 2 , m 3 , m 4 , m 5 ) 
I 
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2mi dmi 



ml 



log 



l c (m 1 ,m 2 ,m 3 ,m 4 ,m 5 ) 

mi + m 2 + m 3 + m 4 
m 3 + m 4 + m 5 



(31) 



(47r) 3 2mi((mi + m 2 ) 2 

In appendix [A] the explicit expressions are given for necessary J- . Equation 
( |29|) can be transformed to a more suitable form by using the symmetry prop- 
erties of the diagram. Obviously, the diagram is invariant in the following 
permutations of masses: 



2 






3 






1 




5 


S)- 




6 


i)- 




4 


I) 



1 

4 

because they are equivalent to rotations of the diagram. 

Using these symmetries, Eq. (EDf) can be written in matrix form 



(32) 



M(mi, m 2 , m 3 , m 4 , m 5 , m 6 )K 



(2ml 
A 
V B 



A 
2m\ 
C 





I A I + J, (33) 
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where / = (1, 1, 1) T and J = (J±, J 2 , J3) 7 ', 



4 2,2 2 

A = m x + m 2 — m 6 

B = m\ + wig — mj 

C = m 2 2 + m 2 3 -mj (34) 



and 



J 1 — J 2 J 2 + ^3 ^3 
J2 = <-^3 — J3 ~\~ J\ — J\ 

J 3 = Jl- 4 + J\- J 2 4 . (35) 

However, since in the following all masses except m\ will be fixed and it 
occurs in the matrix only in quadratic form, the following notation will be 
used for simplicity: 

M(m\) = M(mi,m 2 ,m 3 ,m4,m 5 ,m e ). (36) 



3.2 Degenerate case 

Consider first the degenerate case. If \M(m\)\ = the matrix has an eigen- 
value of zero. Let U be the corresponding eigenvector. Then 



= U l M{mX)K = U 1 II A + U 1 J. 



(37) 



Let Di be the minor of the determinant of matrix Ai(mf) with one of the 
rows and column i removed. Then a suitable eigenvector is 



U 



-D 2 
V ^3 J 



(38) 



Hence, the result for X4 can be written down: 

U T J \Mj\ 



1a = 



U T I 



(39) 



where M.x denotes matrix Ai(m1) with one of the rows or columns replaced 
by vector X. Not every choice of row or column is always possible, but one 
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can use the symmetries of the diagram to transform the matrices to a form 
in which the chosen row is the first one. Therefore the first row will be used 
here. 



As an example let us calculate the special case m% 
777.2. Then the matrix M. reads 



M 



/0 






2m\ 
m\ + 7773 - 



777, 





7772 + m 3 ' 

2ml 



777, 



0, 777 5 

\ 
/ 



777 3 , 777 6 



Then 



Jl 


J2 


h 







2m\ 


m\ + 777,3 — 


ml 





m\ + ml — 777,4 


2m 2 3 




1 


1 


1 







2m 2 2 


m 2 + m 3 — 


m 2 4 





m\ + 7773 — 7774 


2m 2 3 





-Jl 



-J? + J2 



J?. + J3 



•2 ~ °2 °3 

= -X C (m 2 , 777 3 , 777 2 , 777 3 , 777 4 ) + Xc(777 2 , 777 4 , 0, 777 3 , 777 3 ) 
-2c(t77 3 , 777 2 , "73, 777 4 ) + Xc(^3, W 4 , 0, 777 2 , 777 2 ) 

Substituting the expression (|3lD one obtains 



1 A {0, 777 2 , m 3 , 7774, 7773, 7772) 



1 1 



(4tt) 3 2 



log 



2m 2 



"12 +m-j +ni4 



+ 



log 



2m 3 



m2+m3+m4, 



m3(m2 ~ ( m 3 + m D) 7772(7773 — ( m 2 + m D) 



, 1 1 
+ ( — + — 

,777 2 7773 



log 



2(m2+m 3 ) 

7Tt2+m3+7Tt4 



7772 



(ml + 777§) 



In the case 7772 = 7773 = 777, 7774 = 
Z/i(0, 777, 777, 0, 777, 777) 



this result simplifies to 
1 1 , 

-log 2), 



(47r) 3 4?77 3 



(40) 



(41) 



(42) 



(43) 



which is the result of 25 
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3.3 Non-degenerate case 

Consider now the non-degenerate case \ A4\ ^ 0. Now the matrix is invertible. 
The equation can then be solved for K\. This leads to a first order linear 
differential equation: 



dm\ \M(m\)\ 

where in M.x substitution of X to the first row is assumed. From now on 
for Ia a notation similar to (|36|) will be used: 



Z A {m\) = T A (m 1 ,m 2 , m 3 , m 4 , m 5 , m 6 ). (45) 
In the domain where the coefficients on the right hand side are continuous 



(fPj) now has a unique solution. If Ia is this solution, it has a form 



I Aim,) = -exp(-/ -^dt 



6XP( ' \M(t)\ dt ) \M(s)\ ds + C 



(46) 



The value of the integration constant C is to be determined. Let the 
lower limit of the integrations be m 2 , a point such that the coefficients are 
continuous. Now let m\ approach m 2 to see that in this case C = X^(m 2 ), 
since the exponential factor approaches unity and the integral in the brackets 
vanishes. Now let m be such that |*M(mQ)| = 0. Since /^(mf) coincides 
with lA{m 2 ) when m\ < Wq and Ia is continuous also when ml = ?71q, they 
must coincide also at the point m\ = for Ia to be continuous. Thus, 
when m\ < Wq, 

lA(ml) = lA(ml). (47) 

However, when ml — > tbq, the exponential factor in Q4"6] ) diverges. Therefore 
the expression inside the brackets must vanish and one obtains 

i A (m\) = - r' -p f /* W^"- («) 

v v Jmi \Jmi \M(t)\ J \M(s)\ v ' 

Now the following relation is true: 

\M I {ml)\ = \-^\M{m\)\. (49) 
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This can be seen as follows. Consider the determinant expanded as a sum of 
minors: 



\M{m\)\ = 2m\ 



2ml 


C 


-A 


A 


C 


+ B 


A 2m\ 


C 


2m\ 


B 


2m\ 


B C 



(50) 



The first minor is constant in m\ and dA/dm\ = dB/dm\ = 1. Therefore 



d 



dm\ 



\M{m{)\ 



2m\ 



C 



2m\ 



-A 



1 

A 
B 



C 
n 

1 C 
1 

7 

c 



1 2m\ 



+ B 



A C 
B 2m\ 
1 2m\ 



+ B 



A 2m\ 
B C 



1 



A 2m\ 



1 

C 



+ 



A 2m\ 



2m\ 
1 

c 



c 

A 



1 2m\ 



B 

C 



2m\ 



1 

2 

1 C 
2\Mi(ml)\. 



2m\ 



(51) 



1 

B C 

Thus, the integral in the exponential can be calculated and one obtains 

'\Mj(s 



1a« 



\M(s) 

ml \ \M{m\)\ \M{s)\ 

I \Mj{s) 



\M{m\)\ 
2 



-ds 



zds 



\M(s)\ 
-i |^j(x 2 )| 



yJ\M(ml)\ Jm ° ^J\M(x 2 )\ 



xdx. 



(52) 



This result holds only if the integrand has no singularities between m and 
mo. However, Aij diverges only if mi = 0, m2 = or 777,3 = 0, i.e. there is 
no closed massive loop in the diagram. That case must be treated separately 
and leads to infrared divergences. An appropriate choice of mo ensures that 
the denominator Ai(x 2 ) has no zeros on the domain of integration. 

When evaluating the integral (^) one needs the following integral 



log(l + at)dt 



Vl-/9t 2 + 7^ 4 (l + St) ' 



(53) 
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Unfortunately this integral is not expressible in terms of usual special func- 
tions. Therefore one needs numerical evaluation of the last expression. Since 
the integrand consists of only elementary functions, this can easily be done. 

As an example consider the case m\ = m 2 = m 3 = m 4 = m 5 = m 6 = m. 
Now |.M(x 2 )| = 6m 4 x 2 — 2m 2 x 4 and 

\Mj{x 2 )\ = 3m 4 -x 2 m 2 (J 2 + J 3 ) 
m 



3 



4 3m + x 
log log 

3 2m + x 



(4< 

x ( x 4m 3m + x 

+77 — ; — 77 lo S 77 — ; lo S 



2m + x\ 2m — x 2m + x 3m 
Since m can be chosen to be zero, the result is 

X^(m, m, m, m, m, m) = 

1 If 1 dx / 3 , 3+x 

/ — i^=^= lOE h lOE 

(4vr) 3 m 3 v/2 Jo y/S^x^K 4 6 2 + x 
x 2 , 4 x , 3 



(54) 



lo g^— + 77-— log^T- • (55) 



4-x 2 & 2 + x 2 + x 
Numerical evaluation of this integral gives 

T t \ °- 0217376 

XAm, m, m, m, m, m) rs - . - — — . (56) 

(47r) d m d 

3.4 Master integral 

The integral X a is called the master integral, since all the other two loop two 
point integrals can be obtained from it by removing some of the propagators 
from the integrand. Since X a is the two point counterpart of X4 and it 
converges, it can be calculated using Lemma [TJ: 

I a (p;m 1 ,m 2 ,m 3 ,m 4: ,m 5 ) = 
2yri 

— (I A {m 3 , mi, m 2 , ip, m 5 , m 4 ) - X A (m 3 , mi,m 2 , -ip, m 5 , m 4 )) . (57) 
V 

In case of two point diagrams it is more convenient to use the following 
notation 

M(p; m lj m 2j m 3j m 4 ,m 5 ) = M(m 3 ,m l ,m 2 ,ip,m 5 ,m 4: ) 
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/ 2mJ ml + m\ — m\ m\ + m\ — m\ \ 

= m\ + m\ — m\ 2m\ m\ + m\ + p 2 (58) 

^ m 2 . + m| — m\ m\ + wij + p 2 2m 2 , J 

and again, for simplicity 

M(p; ml) = M(p; mi,m 2 , m 3 , m 4 , m 5 ). (59) 

In the matrix Ai only the squares of masses are present and therefore those 
are equal for both terms. One obtains 

rr t 2\ 2 f 113 , 1 

l a {m 3 ) = = / Gte 



|A^O;m|)| ^|.M(p;a; 2 )| 
(|A^j(x,mi,m 2 ,ip, m 5 ,m 4 )| - \M J (x,m 1 ,m 2 , -ip, m 5 ,m 4 )|) . (60) 

Now the matrices A4j differ only in the first row. Hence 

— (\Mj(. . . ,ip, . . 01 - |-Mj(- ■ ■ , -ip, ■ ■ 01) = \M H {. ■ ■ ,ip, ■ ■ 01, (61) 
p 

where 

27T1 

//(.... i/,...) _(./(....!/,...) ./(.... -i/,...)). (62) 
p 

Since vector J consists of convergent vacuum integrals, one can use Lemma 
l| to notice that H consists of corresponding two point integrals: 

H 1 (p;m u ...,m 5 ) = Hf - R\ + E\ - E\ 

H 2 (p;m u ...,m 5 ) = + - H* 

F 3 (p;m 1 ,...,m 5 ) = tf 3 2 - tff + 2 , (63) 

where Hf is, similarly to J- , integral If or Ih which has been constructed by 
removing the particle j from integral X a and taking square of the propagator 
of particle i. The explicit expressions of these functions are given in appendix 

B 

If one writes 

Mh{p; mi, m 2 , m 3 , m 4 , m 5 ) = Mii(m 3 , mi, m 2 , ip, m 5 , m 4 ), (64) 
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Eq. ( pO| ) can be written in the form 

I a (p; mi, m 2 , m 3 , m 4 , m 5 ) = 

2 r m 3 \J\AH(p',rn 1 ,m 2 ,x,m 4: ,m 5 )\ 



f m3 \jyiH{P;m 1 ,m 2 ,x,m4,m 5 )\ ^ 



\M(p; ml)\ Jm o ^J\M(p;x 2 

where m is again a root of the equation 

\M(p;m 2 )\ = 0. (66) 

As an example, consider now the case m\ = ■ ■ ■ = m^ = m. Then the 
matrix A4 (p; x 2 ) = M. (p; m, m, x, m, m) is 



M(p; x 2 ) 



( 2x 2 x 2 x 2 

x 2 2m 2 2m 2 + p 2 \ , (67) 
\ x 2 2m 2 + p 2 2m 2 



and the functions f/j are 



(Air) 2 mp 2 (p 2 + 4m 2 ) 



( p p \ p +(2m + x) 
p arctan arctan — m log — — — , . 

V 2m 2m + xJ (2m + x) z J 

1 / p p \ 
H 2 = H 3 — — r arctan arctan I . (68) 



(AiT) 2 2mpx(2m + x) \ 2m 2m + x 

Substituting these to Eq. (EST) gives the result 



I a {p'i m i m i m i m i m ) — 

1 r m dx 



(4tt) 2 mp 2 ^p 2 + 3m 2 Jo J( p 2 + 4m 2) _ % 2 



2p ( p p \ p + (2m + x 

arctan arctan + log — — 

2m + x \ 2m + x 2m J (2m + x' 2 



,2 n 



(69) 



4 Effective potential and self energy in scalar 
theory 
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Figure 5: Feynman rules of the scalar theory 

4.1 Lagrangian 

As an application of the integrals calculated in the previous section scalar 
A</> 4 theory will now be discussed. This is the simplest possible nontrivial 
quantum field theory. However, it has also physical significance. In theory 
of critical phenomena it is in the same universality class as the Ising model. 
The Lagrangian of the theory is 

C = \{d4? + \mltf + \\<t ) \ (70) 

Suppose now that m\ < 0. Then the minimum of the Lagrangian is not 
anymore in the origin. Now make a shift — > 0o + <P to get a new broken 
Lagrangian. The terms linear in fields are discarded, since they cancel the 
tadpole terms in the true minimum when calculating the self energy. When 
calculating the effective potential they are also discarded, so that a presen- 
tation in terms of vacuum diagrams can be obtained |?TJ. The Lagrangian 
now reads 

£ = \™l<t>l + \tyo + \^<Pf + \™ 2 4> 2 + A0o0 3 + \\<P\ (71) 

where 

m 2 = m 2 + 3\<p 2 . (72) 
The Feynman rules of the theory are shown in Fig. [5| 

4.2 Self energy and renormalization 

Let us start analyzing the system by calculating the self energy. In Fig. |6] 
it is expanded to two loop order. The number of different Wick contractions 
corresponding to each diagram is given in table [I]. 
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a 


1 


a 


15552 


A 


1296 


P 


36 


b 


15552 


B 


1944 


7 


6 


c 


864 


C 


216 






d 


432 


D 


216 






e 


288 


E 


24 






f 

g 


864 
192 


F 


72 






h 


432 







Table 1: Contraction numbers of diagrams 



In dimensional regularization only diagram g diverges. To remove this 
divergence a mass counterterm 5m 2 must be introduced. Using the given 
symmetry factors and integrals, the value of the divergent diagram is 

6A 2 / 1 3 3m p 1, fi 2 \ 

' 77 + 7 ~— arctan 777 + 7 lQ g n^ 2 , 77 ■ (73) 



(4n) 2 \4e 2 p 3m 2 to 9m 2 + p 1 

The correct value of the mass counterterm is then 

*- 2 - (54- (74) 

This gives rise to a running mass 

m ^ ) = (Sp log £:' (75) 

where A m is a dimensional parameter such that mg(A m ) = 0. Thus, the 

mass m 2 is a function of both the renormalization point and the vacuum 
expectation value of the field: 

m 2 = m 2 (7i,0o)- (76) 

Since no other divergences are present, the coupling constant A does not run. 

The self energy can now be calculated by collecting all the relevant inte- 
grals and the two loop part of the result is 
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(4tt) 2 I 

+A 3 2 

54 
pm 



A 



9 m p a 

- - 18— arctan — + 3 log — — — ^ 

2 p 3m 9m z + p z 

54 



9 54 / p \ 

f QjJ^ I 

p 2 + 4m 2 m 2 k 2 \ 2m J 



2 log 3 arctan h i I Li 2 

2m 



V I ( l P \ ( 2m — ip 
1 + Li ' 



3m J 



L12 f 7t- ) ~ Li 2 



\3mJ 



2m + ip 



in 



27 



., , 1/'('V • Llm 2 ) arctan 

my(p z + 4m' : ) V 3m 

P 



+ (6 log 3 — 8 Wp 2 + 4m 2 ) arctan — 6m(p 2 + 2m 2 ) log 1 + 

2m y 

P-u 2 (p 

3mJ \3mJ 

2+*)-U,(-2-* 

my V m 



9m 2 



+3ip(jo 2 + 4m 2 ) ( Li 2 
Li 2 
4 



648 



dx 



mp 2 yjp 2 + 3m 2 Jo ^(p2 + 4m 2) _ x 2 

( 2p ( p 

arctan 

y2m+x V 



p \ p 2 + (2m+x) 

arctan + log — 7 — — 

2m+x 2m/ (2m+x) 2 



.(77) 



4.3 Effective potential 

The calculation of the effective potential is very similar to that of the self 
energy. The diagrams to three loop order are given in Fig. [7| and the numbers 
of contractions in table |l| 

The result can be written down at once 



m 



+ 



+ 



(4tt) 2 
1 



Am 2 - 3A 2 0o 



(4 



7T 



4 
mA 2 



1 A* 1 
log — + - 

5 3m 2 



P 27 

31og 4^ + ¥ 
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n(A) = 



'24 



1 + l 

'2 








Figure 6: Loop expansion of the scalar self energy 



m 
A^f 
m 3 



9 9 5 

— + -7T Z 

2 4 



27 



log- 



27Li 2 



27 2 81 / 4\ 2 cj(1 4 81 /1\ 
- y ^ + T (log-) +541og- + yLl2 y 



-54 



3 , 3 + x 
log - + log — — 

4 2 + x 



x 4 a; 3 + x 
log o 1" o io S — n — 



X' 



X 



X 



(78) 



This result agrees perfectly with that evaluated numerically in |26|. 



5 Conclusions 



In this paper, the integrals necessary for the self energy to two loop level 
in a three dimensional scalar field theory have been evaluated explicitly as 
well as the ones necessary for the effective potential to three loop level. In 
almost every case the result can be expressed with elementary functions and 
dilogarithms. 

A large part of the paper has been devoted to the evaluation of the two 
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V(cf> ) = 




Figure 7: Loop expansion of the effective potential 



most difficult integrals, the master integral and the Mercedes integral. Those 
are expressed in terms of a one dimensional integral representation with an 
integrand consisting only of elementary functions. This form makes numeri- 
cal evaluation easy. 

The results have been applied to A0 4 scalar theory. It will be very inter- 
esting to extend them to gauge theories with scalars, like the U(l)+Higgs 
or SU(2)+Higgs models. This will lead to a large number of new diagrams. 
However, all the two loop integrals contributing to the self energy in a gauge 
field theory can be decomposed to a sum of scalar integrals. This is a labo- 
rious task for which a computer algebra system is needed. The results will 
help to deepen our understanding of the phase transitions in gauge theories, 
for example the electroweak phase transition of the early universe. 
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A Integrals 

This is a complete list of the integrals corresponding to the diagrams of Figs. 
H 1 and |. 



where 



and 



at) 

0) 
7) 
a) 



m 

47T 



I ; : ( 2 + 2 log ^- 



Airp 
1 



arctan • 



V 



mi + m 2 
— + - + log 



(An) 2 \4e 2 



mi + m 2 + m 3 



^\M(p;m 1 ,m 2 ,m 3 ,m4,m 5 )\ 
m3 \ Mh(p\ mi,m 2 , x, m 4 , m 5 ) \ 
"'" J\M(p; mi, m 2 , x, m±, m§)\ 



xdx, 



M(p; mi, m 2 , m 3 , m 4 , m 5 ) 



/ 



2m§ 
\ m\ + m| - 



m{ + mj — m\ m\ + wij — m\ 



m\ 
2 

mi 



2m\ 
m\ + m\ + p 2 



m\ + wij + p 2 
2m\ 



(79) 
(80) 
(81) 



(82) 



(83) 



M H {p; mi, m 2 , m 3 , m 4 , m 5 ) = 
/ #1 H 2 
m\ + m 2 — m\ 2m\ 



\ m\ + m\ 



m 



m\ + wij + p 



The functions mi, m 2 , m 3 , m 4 , m 5 ) are 



H 3 

m\ + m| + p 2 
2ml 



Hi(p;mi, 
H 2 (p; mi, 
H 3 (p; mi, 



,m 5 ) = Hf-Ht + Hl-Hl 
,m 5 ) = H l 2 +Hl-Hl 
m 5 ) = Hl-Hl + Hf, 



(84) 



(85) 



where 



A? 



H x (m 1 ,m 2 ,m 4 ,m 5 ) 
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= H Y (m 1 ,m 3 


m 5 ,m 2 ) 


H\ 


= H x (m 2 ,m 1 


m 5 ,m 4 ) 


Hi 


= H Y (m 2 ,m 3 


m 4 , mi) 


H\ 


= H z (m 4 ,m 2 , 


m 3 ,m 5 ) 


Hi 


= H z (m 4 ,m 3 , 


m 2 ,m 5 ) 


Hi 


= H z (m 1 ,m 3 , 


m 5 ,m 2 ) 


Hi 


= H z (m 5 ,m 3 , 


nix, 777.4) 


Hi 


= H z (m 2 ,m 3 ,m 4 ,mi) 


Hi 


= H z (m 5 ,m 1 , 


m 3 ,m 4 ), 



H x {m 1 ,m 2 ,m 3 ,m 4 ) = 

1 1 p 

■ arctan ■ 



(47r) 2 2mip ((mi + m 2 ) 2 + p 2 ) m 3 + m 4 

H Y (m 1 ,m 2 ,m 3 ,m 4 ) = 

1 1 
(47r) 2 2m\p ((p 2 + m 2 + m 4 ) 2 — 4m 2 m 2 ) 

(^(p 2 + m\ — m\) arctan 



+777 ip log 



m 1 +m 2 + 7773 

P 2 + (mi + 777 2 + m 3 ) 2 \ 



(777 2 + 777 3 + 777 4 ) 2 J 

H z (m 1 ,m 2 ,m 3 ,m 4 ) = 
1 1 

(47r) 2 2m 2 p{m\ — [m 2 + m 3 ) 2 ) 

( p p 
arctan arctan ■ 



V 777! + 7772 + ?™3 m l + m 4 



f , 7772+7773+7774 

<^ Z777i log arctan 

I 777 2 + 777 3 — 777 4 



(47r) 2 4p777i777 4 (777 4 — 777 2 ) I 777 2 + 777 3 — 777 4 

777! + 777 2 + 7T7 3 p 

+2m 4 log arctan ■ 



m 2 + m 3 — 777! + 777 5 

\ T . ( 777 5 + 777i - i]9 \ T . / 777 5 + 777i + ip 

+17T7 4 L12 ; + Ll 2 ; 

V 777 2 + 777 3 — mi) \ m 2 + 777 3 — 777i 
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c) 

d) 

e) 
/) 



9) = 



h) 
A) 



- Li 2 
+imi 

- Li 2 



m 5 - mi - ip 



Li 



m 1 +m 2 + m 3 

m$ — m 4 + ip 



— Lio 



m 2 + m 3 + m 4 
r«5 + m 4 + ip 



+ Li 2 



m 5 — m x + ip s 

mi+m 2 + m 3l 
m 5 — m 4 - 



\p 



m 2 + m 3 — m 4 



m 2 



(47r) 2 p(m 2 — m§) V 

(47r) 2 (m 2 -m 2 ) bg 
1 m 2 



( arctan 



- Li 2 

P 



mi + m 4 
^2 + m 3 + m 4 
mi + m 2 + m 3 ' 



^2 + m 3 + m 4 
ms + m 4 — ip 

"^2 + — TTl4 
p 

— arctan ■ 



m 3 + m 4 



{An) 2 mi + m 3 



1 



(47r) 2 4pm 4 



2 log 



m 2 + m 3 + m 4 



+i 



Li 5 



+ Li 2 
1 

(I^j 2 
1 



m 2 + m 3 
mi + m 4 — ip 
m 2 + m 3 — m 4 
mi — m 4 + ip 



arctan • 



p 



m 4 
- Li 2 



— Li? 

m2 + Tn 3 + 
1 3 m x + m 2 + m 3 
45 + " 



m x + m 4 
mi — m 4 — ip 
m 2 + m 3 + m 4 
mi + m 4 + ip 



m 2 + m 3 
arctan • 



- m 4 
p 



P 



mi + m 2 + m 3 



log 



fi 2 



(mi + m 2 + m 3 ) 2 + p 2 / 
p 



(47r) 2 p : 



arctan • 



mi + m 2 
2 



arctan • 



p 



m 3 + m 4 



^\M(mi, m 2 , m 3 , m 4 , m 5 , m 6 ) | 
mi |7Wj(a;, m 2 , m 3 , m 4 , ms, mg)| 



mix 
J mo 



^\M{x, m 2 , m 3 , m 4 , m 5 , m 6 ) 
where the matrices .M and .Mj are 



■ cc doc ■ 



7W(mi, m 2 , m 3 , m 4 , m 5 , m 6 
/ 2m 2 



m\ + wij 



m\ + m 2 - 
2m 2 



m 6 mf + m 



mi 



\ m\ + m\ — m\ m 2 , + m\ — m\ 



m 2 + m 3 — m 4 

2m 2 
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Mj(m 1 ,m 2 ,m 3 ,m4,m 5 ,m e ) = 

(Ji J 2 J 3 

m\ + m\ — m\ 2m\ m\ + m\ — m\ 

m\ + ml — ml m\ + m\ — m\ 2m\ 

The functions Jj(mi, m 2 , m 3 , 7714, m 5 , m 6 ) are 

J\ = J2 ~ ^2 J3 ~~ J3 
J2 — J3 — J3 + J\ — J\ 

J3 — J\ — J\ + j 2 ~ J21 

where 



A 


= ic 


[m 2 


m 3 


m 5 


m 6 ,m A ) 


4 


= ic 


[m 2 


m 4 


m 1 


7775,7773) 


A 


= ic 


[m 3 


m 2 


m 5 


m 6 ,m A ) 


7 5 

J 3 


= ic 


>3 


m 4 


m 1 


m 6 ,m 2 ) 


T 2 

J 3 


= ic 


>3 


m 4 


m 1 


m 6 ,m 5 ) 


4 


= ic 


[m 3 


m 5 


m 2 


m 6 ,mi) 


Jl 


= ic 


[mi 


m 6 


m 3 


7774,7775) 


Jl 


= ic 


[mi 


m 5 


m 2 


7774,7773) 


Jl 


= ic 


[mi 


m 2 


m A 


m 5 ,m 6 ) 


7 5 


= ic 


[mi 


m 6 


m 3 


m A ,m 2 ) 


7 3 

J 2 


= ic 


[m 2 


mi 


m A 


m 5 ,m 6 ) 


J 2 


= ic 


[m 2 


m 6 


m 3 


777.4, mi) 



and Xc is the derivative of Ic'- 

i c (m 1 ,m 2 ,m 3 ,m 4 ,m 5 ) = 

1 777i + m 2 + m 3 + m 4 

log 

(47r) 3 2mi((mi + m 2 ) 2 — 7775) 7773 + 7774 + 7775 



n\ 

(47r) 3 4?775777 6 (7776 — 7T7§) 
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C) = 



+2Li 2 
-m 5 

+2Li 2 
1 



log 



m 3 + m 4 + m 5 l (m 3 + m 4 ) 2 — m 2 



log 



rri3 + f^A — "^5 
mi + m 2 - m 5 

m 3 + m 4 + m 5 

m 3 + m 4 + m 6 

log ■ log 

m 3 + m 4 — m6 

mi + m 2 — m 6 



(m 4 + m 2 + m 5 ) 2 
T . ( m 1 +m 2 + m 5 

Z Lily 



(47r) 3 2m 5 



m 3 



log 



- m 4 
m 3 -\ 



h m 6 
m 4 - 



m 3 + m 4 — m 5 
(m 3 + m 4 ) 2 - m| 
(m 4 + m 2 + m 6 ) 2 
T . / mi + m 2 + m 6 



m 5 



m 3 



+Li 2 



mi + m 2 



• m 4 
m 5 



m 3 + m 4 + m 5 



-m 5 
- Li 2 



log 



m 3 + m 4 — mg / J 
(m 3 + m 4 ) 2 — ml) 5 

m x + m 2 + m 5 
m x + m 2 + m 5 



m 3 + m 4 — m 5 



D) 
E) 
F) 



m 5 



(47r) 3 (m| - ml) 
1 4 



log 



mi + m 2 + m 4 
mi + m 2 + m 3 ' 



1 mim 4 
(47r) 3 m 2 + m 3 



.(101) 



(102) 
(103) 

(104) 

(105) 
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Erratum to "Feynman diagrams to three loops 
in three-dimensional field theory" 
[Nucl. Phys. B 480 (1997) 729] * 

Arttu K. Rajantie 1 

Department of Physics, P.O. Box 9, FIN-00014 University of Helsinki, Finland 

(N 

■ The sign of the constant term —2 was incorrect in Eqs. (21) and (A. 26): 

<n : 

! 1 A / 1 1, 2rrii , X^w/ 

^E m > "J- " 2 + olog— 1 + log^^ 
br) 3 ^ V 4£ 2 /i 



43 



X 



(4tt) 3 ^ 'V 4 £ 2 £ /2 

\ One m 2 should be m 4 in Eq. (A. 9): 

O 

43' Hy(mi,m2,m 3 ,m4 



p ' ' ' (An) 2 2mip ((p 2 + m 2 + m^) 2 — Am 2 m 2 ) 

x ( (p 2 + m 2 — m 2 ) arctan 



mi + m 2 + m 3 

I mip log + ( m i + m 2 + ^) 2 \ / 2 x 
(m 2 + m 3 + m 4 ) 2 / 



Five signs were incorrect in Eq. (A. 10): 



■ , m^ma + ms p 
(0) = — 2 2T 1 log arctan ■ 



(47r) 2 4pmim 4 (m 2 — m 2 ) I m 2 + m 3 — mi mi + m$ 

m 2 + m 3 + m 4 p 
— 2mi log arctan ■ 



m 2 + m 3 — m 4 m 4 + m 5 

+im 4 



r . . m 5 + mi - ip \ / m 5 + mi + ip 

L12 ; - Li 2 



m 2 + m 3 — mi / V m 2 + m 3 — m x 
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Li 2 



—mil 



m 5 — mi — lp 



Li 



mi + m 2 + m 3 

m 5 + m 4 — ip 



+ Li 2 



Li 2 



m 2 + m 3 - 
m 5 — m 4 — ip 

m 2 + m 3 + m 4 



-Li 



m 4 
+ Li 2 



m 5 — mi + ip 
m! + m 2 + m 3 

m 5 + m 4 + ip 
m 2 + m 3 — m 4 
m 5 — m 4 + ip 

"^2 + + m 4 



(3) 



In three of the quantities Jj listed in Eq. (A. 21), one or two mass indices were 
incorrect: 



J^=X c {m 3 ,m 1 ,m A ,m 6 ,m 5 ), 
J\ =lc(rni,m 3 ,m4,m 6 ,m 5 ), 
J\ =l c (m 2 , rriQ, m 3 , m 5 , mi). 



(4) 



In Eq. (A. 23) the 7r 2 /12 terms were missing: 



{B) = 



(47r) 3 (ml - ml) 



m 5 



7r 2 If m 5 + mi + m 2 

1 — log 

12 4 V m 5 + m 3 



1 / Ms — mi — m 2 

+ -Ll 2 ; ; 

2 \m 5 



m 6 



m 3 + m 4 

mg + mi 



+ m 4 

1 / Ms — m 3 — m 4 
- Li 2 

2 \m 5 + mi + m 2 



m 2 



t 1 /, 
h - log 

12 4\ mQ + m 3 + m4J 



1 t . fm 6 - mi - m 2 
+ - Li 2 

2 Vm 6 + m 3 + m 4 



1 t . fm 6 -m 3 -m 4 
H — Li 2 

2 \m 6 + mi + m 2 



(5) 



On p. 740 the sign of the integration constant C should be the opposite. On the 
fourth line of p. 742 there should be ^J\M(x 2 )\ and on p. 743, below Eq. (59), 
\M(p;ml)\=0. 

These errors change neither the results in Eqs. (77) and (78) nor any of the 
conclusions. 

I would like to thank F. Eberlein, S. Larin and M. Stroesser for pointing out 
these mistakes for me. 
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